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DEFLECTIONS GRIDWORKS AND SLABS 


WALTER JUN. ASCE, SHIGEO OKUBO,? 
AND JOEL JUN. ASCE 


Because the exact determination elastic deformations plates and 
slabs the mathematical theory elasticity has proved impractical 
all but the simple cases, engineers have been forced devise approximate 
means solution. horizontal gridwork beams intersecting right angles, 
with deformation characteristics analogous those slab, will develop 
elastic surface under normal load similar the surface the slab under the 
same load. technique presented determine the elastically deformed 
surface plate slab that has been transformed into equivalent grid 
system. method employs auxiliary force system for controlling vertical 
displacements the joints and moment and torque distribution process for 
transmission the displacement effects. 


INTRODUCTION 


Approximate methods determining the deflected surfaces slabs are 
not new. system beams, intersecting right angles one another, 
forms gridwork that yields deflected surface similar that slab when 
analyzed under normal loads. This statement appears obvious when one con- 
siders the similar elastic surfaces developed slab and simple window 
screen supported and loaded like fashion. The analysis such gridwork 
not simple, however. 

One method determining the deflection gridwork beams pre- 
sented Miklés who assumes that the individual beams comprising 
the gridwork deflect without rotation their intersections with other beams. 

comments are invited for publication; the last discussion should submitted 
March 1952. 
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Simultaneous differential equations are then solved order obtain joint 
displacements. 

When assumed that bending twisting moment transmitted 
grid beam intersections, another method grid analysis possible. With 
such assumption, Stephen Timoshenko employs trigonometric series 
express the curves developed the individual grid grid- 
work, unable transmit bending and torsional moments, not, however, 
analogous slab plate its action. 

realistic rigid joint, capable resisting bending and torsion, assumed 
each beam intersection, the analysis the resulting grid tedious process. 
Linear equations can formed involving three unknown elements each 
joint—the loads the beam, the moment change the longitudinal beam, 
and the moment change the transverse Both moment changes are 
developed torsional effects the mutually perpendicular beams. 
inordinate amount time involved the solution the simultaneous 
equations particularly since there are three times many equations there 
are joints the grid. 

technique here presented for eliminating two the three unknown 
quantities each joint moment and torque distribution process. 
some its aspects, this procedure similar the method finite differences 
developed Henri Marcus’ and extended skew slabs Vernon 

Finite difference equations have been applied rectangular and skew 
slabs divided into network points determined arbitrarily ruling the sur- 
face with series intersecting lines. Moments and deflections are obtained 
each point the network solving simultaneously two groups linear 
equations: one group containing moment values unknowns, the second group 
expressing deflection values unknowns. The number equations 
solved one group is, the general case, equal the number points the 
arbitrary network. The solution the two groups simultaneous equations 
can avoided the original fourth order differential equation the plate 
not divided into two parts. When finite difference equivalents are determined 
directly from the fourth order equation, solution required for only one 
series linear equations. Nevertheless, this procedure not 
Mr. Marcus. 

comparison, the method used this paper applicable physically 
analogous grid rather than the ruled surface the slab itself. After the 
bending moments and torsional moments have been distributed over the grid, 
only one series linear equations need solved order define the deflection 
pattern. This possible because the equations are written terms un- 
known deflections produced auxiliary forces the grid points. Bending 
moments and torsional moments can then found without recourse the 


die Biegung von Stephen Timoshenko, Zeitschrift fiir angewandte Mathe- 
matik und Mechanik, Band 13, 1933, 153. 

und Platte,’’ Béla Enyedi, Beton und Eisen, Vol. 34, Heft 1935, 42. 

Theorie elasticher Gewebe und ihre Anwendung auf die Berechnung biegsamer 
Henri Marcus, Julius Springer, Berlin, 1924. 


Skew Slabs,’’ Vernon Jensen, Bulletin Series No. 332, Univ. Illinois, Urbana, 
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second series simultaneous equations that are required the simplified 
finite difference solutions. 


PROCEDURE 


The procedure for determining the elastic surface plate slab sub- 
jected normal loading consists the following steps: 


The arbitrary division the slab into orthogonal strips that are con- 
sidered the beams the analogous gridwork. 

The determination factors, based the individual beam’s resistances 
bending and twist, for distributing the unbalanced moments each 
grid beam joint. 

The vertical displacement each singular joint the grid and conse- 
quent introduction fixed-end moments the beams. Joints where deflec- 
tion produces distinctive moment and torque pattern are known singular 
joints. 

The conversion the fixed-end moments into bending and twisting 
moments moment and torque distribution process. 

The computation reactions all joints terms the unknown dis- 
placements. 

The computation the moments and torques each joint which were 
expressed the distribution process terms the arbitrary displacements. 


Transformation the Slab into Comparable plate slab, when 
divided into strips two mutually perpendicular directions, can trans- 
formed into grid system provided that each strip considered grid 
beam. fact, analogy between slab and gridwork indicated the close 
correspondence slab deflection values found the grid method and checked 
the theory elasticity. Obviously, finer division the slab will increase 
the number grid beams and will lead more accurate definition the 
elastic surface. 

the slab shown Fig. divided into two series parallel orthogonal 
strips, the resulting system simulates group rigidly connected continuous 
beams mutually supported along their centerlines. For example, the east- 
west and north-south grid strips are considered continuous over supports 
1-1, 2-2, 3-3, and over supports located the centerlines 
the two perpendicular strips. Furthermore, beam segment such 
assumed free rotate and twist over its supports. Along the edges 1-2 
the resistances bending and twist that are presumably neglected the 
east-west beams are actually considered vertical shear forces the north- 
south beams. Because the cross section any grid beam identical that 
comparable slab strip, the gridwork actually pseudogrid. The over-all 
general appearance the system different from that the slab because 
openings exist between grid beams. 

Since the equivalent gridwork loaded and consequently subjected 
moments and torques, relative stiffnesses and distribution factors the 
individual grid beams must calculated. 


\ 
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Relative Stiffnesses and Distribution applied moment (M), 
the plane the beam released the rigidly connected system shown 
Fig. the tangent the elastic curve will rotate through angle 
torsional fixity realized points and the beam will twist through 
the same angle point Similarly, moment applied the plane 
will produce rotation the tangent the elastic curve that plane and 
consequent twist the beam 


Center Lines 


order resist the applied moment bending moments and Mic 
will developed beam and torsional moments and will 

When the rotation through the angle has taken place, 


For unit rotation beam using slope deflection equations, 


For unit twist the beam with circular cross section: 


GI, 


with rectangular cross section, or, 
3 


ticity shear; polar moment inertia beam grid strip; moment 
inertia with respect axis perpendicular the plane bending; 


Nostrand Co., Inc., New York, Y., 1949, pp. 265-266. 
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elasticity. With the expressions for moment the longitudinal and transverse 
beams known, distribution factors can written for any joint. For example: 


Distribution factors for the ordinary grid will the same every point 
when lengths between intersections are made constant each direction and 
identical cross sections are selected. 

Force Systems for Controlling moment distribu- 
tion process used cannot applied directly any loaded structure whose 
joints undergo translation. order control prevent these joint transla- 
tions, necessary calculate series auxiliary force systems that arise 
from the arbitrary displacements joints. These forces are applied every 


= 
a 4 


DEFLECTION JOINT 


joint, and the translations they cause prevent must such produce 
equilibrium with the forces actually acting the structure. 

After the distribution factors for the unbalanced moments and torques 
every joint have been determined, arbitrary displacement (A) introduced, 
turn, each singular joint shown Fig. singular joints the 
particular grid illustrated are points and all other 
joints will induce similar moment and torque patterns. The auxiliary force 
system shown consists reactions all joints other than joint which 
point the displacement has been produced. These reactions prevent dis- 
placement from the horizontal datum plane but not prohibit rotation and 
twist the grid members. The lengths the grid members being equal 
both directions, the deflection induces fixed-end moments equal 
the initially fixed joints and joint where the deflection occurs. 

Moments are placed the grid first the east-west direction and are 
assigned the customary moment distribution signs from the south. 
Ordinary fixed-end moments need only introduced the beams one 
direction, for the symmetry the gridwork yields identical distribution pat- 
terns the mutually perpendicular directions. 


| 
8 
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Next, the joints are progressively released and balanced, introducing bend- 
moments and torques all beam intersections. After all the joints are 
balanced, carry-over process effected, torque moments carrying over 
opposite ends the beam their full values but with opposite sign, and 
bending moments carrying over half values with the same sign. The second 
cycle completed rebalancing each joint with the usual distribution factors. 
Three four cycles generally furnish thorough dispersion the original 
fixed-end moments. 

When the distribution process complete any one joint, bending moment 
values the beam ends are used determine the reactions terms the 
displacement. Although distribution performed the direc- 
tion only, similar group bending moments must placed the 
south beam order define the complete reaction pattern. Fig. shows 
assumed deflection pattern due displacement joint The notation 
should self evident. For instance, the reaction joint caused 

deflection joint should remembered this point that all 
reactions are expressed coefficients the particular grid 
selection, similar displacements must introduced the other singular points 
and and reactions must calculated. 

Using relationship noted that the reaction any point 
for example, due the deflection will equal the reaction 
due the displacement joint 


Thus the complete reaction joint with all points the grid displaced 
through unknown values, will be: 


Since only singular points exist, joints and 


lar equations can written points the grid and equated load that 
might exist the intersection, else equated zero there load 
the joint. 

the case uniform load across the plate, load exists every inter- 
section and the displacement three points will define the elastic surface. 
Consequently, this simple case, only three simultaneous equations need 
solved fur the unknowns and 

With the deflection values known each point the moments and torques 
expressed the distribution process terms can accumulated 
simple multiplication and addition. 
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Example A—Deflections Simple this first application 
the analysis attempt has been made draw analogy between gridwork 
and slab. Fig. illustrates gridwork consisting circular bars rigidly 
connected each junction, with complete fixity 
against bending and twist realized beam ends. 
will determined this structure when 
subjected equal normal loads all joints 
and when single load applied joint 
Since all lengths are equal and cross sections are 
the same, the distribution factors each the joints 
are identical. The distribution factors are com- 


Assuming joints and temporarily fixed 


0.72 0.72 


which unbalanced moment released the joint. Likewise, 


Similar equations can set for the remaining points the gridwork. 
The percentage the unbalanced moment resisted bending moments 
developed beams 42.37%, and the percentage resisted torque moments 
the north-south beams 7.63%. 

are induced points and along the beam and also and 
along the beam shows the fixed-end moments multiplied 1,000 
but exclusive the constant quantity along the east-west beams, 
together with the ensuing moment and torque distribution. The signs the 
moments are accordance with the moment distribution convention when 
viewed from the south. Moments are not distributed the north-south 
direction since they are identical those developed the east-west beams. 

Equilibrating reactions are derived from the end moments the east-west 
direction obtained the The equations equi- 
librium are used determine these reactions after the beams are subdivided 
into their component lengths, shown Fig. Because the constant 
beam length the reaction values shown are multiples the term 
When equal series reactions derived from the north-south end moments 
are superimposed upon the grid and algebraically added the first group 
reactive forces, the equilibrated system shown Fig. results. 
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the theorem reciprocity, the reaction 13.188 junction 
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similar reasoning employed all 


grid points, the total reaction joint can written follows: 
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every grid joint supports load the following expression can ob- 
tained equating vertical forces zero: 


But with equal forces located all joints the symmetrical grid: 
Therefore: 
transposing 


Fig. illustrates the deflected grid subjected equal loads all joints. 


3 
The deflections are expressed coefficients 


When single concentrated load applied joint the deflections 
developed joints and are not equal. equating vertical forces 
all grid joints zero the following linear equations are developed: 


Solved simultaneously, the equations yield the following values for deflection: 


Fig. 8(c) shows the deflections coefficients each the four joints. 
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Since the bending and the torsional moments have been determined 
terms the various displacements, simple matter calculate the mo- 
ments after the deflections have been determined. 

From Fig. the moment due displacement joint 
5,715 and the torque along beam 103 units, 
shown The torque 103 could calculated directly dis- 
tributing fixed-end moments the north-south direction, but this procedure 
unnecessary due tosymmetry. Furthermore, with displacement joint 


(a) ASSUMED ELASTIC SURFACE DEFLECTION SURFACE (c) DEFLECTION SURFACE 
AFTER DISTRIBUTION LOAD JOINTS 1,2,3,4 LOAD PAT JOINT 
MOMENTS DUE 


the moment 1,335 shown point the figure would develop 
1,335 point The change sign arises because original 
fixed-end moments should introduced 6,000 With displace- 
ments all joints the total moment is: 


Substituting values for previously determined, 


which the positive sign based moment distribution convention. 
other moments and torques can determined like fashion. 

B—Analysis Slab, With and Without Center 
second illustration the applicability the method, the 1/16 in.-thick slab 
shown Fig. has been selected for analysis. combination fixed and 
free boundaries has been chosen for the slab order prove that complexities 
arising from variable edge conditions are easily handled. The effects con- 
tinuity will studied determining slab deflections both with and without 
center vertical support joint 12. 

the slab shown divided into series 3-in. wide strips, each strip can 
considered beam equivalent grid system. The grid beams and their 
centerlines, along which the deflections will determined, are shown Fig. 
Because the fixed- and free-end conditions, three sets distribution factors 
are required for the distribution moments and torques over the gridwork. 

Interior joints require but one set distribution constants, for when 
displacement introduced two equal length and equal cross-section beams 


0.075 
0.131 4.015 
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twist through angle two identical but perpendicular beams rotate through 
the same angle. The same stiffness factors are developed, and consequently 
the same distributive effects are required for the distribution fixed-end 
moments the perpendicular direction. Edge joints require different con- 
stants the two directions, for when moments are distributed int north- 
south direction, one east-west beam subjected torque; and when moments 
are distributed the east-west direction two north-south beams are subjected 
twist. The values the distribution factors for interior and exterior joints 
are shown constants can computed assuming all beams 
constant length and also assuming that 0.377 and that the 
ratio beam depth beam width 1/16 in. in. 1/48. The expres- 
sion for moment required produce unit angle twist rectangular 
beam has been taken from Mr. 


INTERIOR JOINT EXTERIOR JOINT 


PLATE 


When the slab subjected uniform vertical load per in., 
the equivalent gridwork loaded each joint with concentration 
Twelve deflections, through the numbered joints shown Fig. 
will completely define the elastic surface. Consequently, equal but arbitrary 
displacements are introduced each the similarly numbered singular grid 
joints turn, while the remaining joints are temporarily fixed horizontal 
datum plane vertical reactions. When all similarly numbered joints are 
displaced simultaneously, symmetrical distribution results and the 
tion can confined one quarter the slab. 

Inasmuch the numerical process distribution similar all joints, 
the procedure illustrated for displacements joint only, shown the 
chart This chart, showing only the north-west quarter the slab, 
facilitates the simultaneous distribution the north-south 
taneous distribution necessitated the different end conditions the two 
directions. For convenience, moments are represented the constant coeffi- 
cient multiplied 1,000. example the carry-over factors for bend- 
ing moment and torsional moment are 1/2 and —1, respectively. The singular 
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joint reactions, exclusive caused the displacements are 
illustrated Fig. 12. 

Similar distributions are performed each the singular joints, each 
operation reducing two the number unknowns. The complete process 
reduces from the number linear equations solved simultane- 
ously. 


After the effects the displacements each the singular joints have 
been determined, the reactions shown Table are calculated, using the 
methods statics. The reactions caused are illustrated Fig. and 
entered Table 1(a). 

When the reactions—expressed terms the unknown displacements 
but still exclusive the term accumulated each joint and each 
reaction equated vertical load 1,000 the twelve equations used deter- 
mine the twelve A’s result. Since the equations are merely restatement 
Table they are not repeated this point. 

With center support joint 12, only the first eleven equations are re- 
quired, inasmuch known equal zero. Furthermore, the 
term not present the eleven equations solved. 

The solution these equations facilitated using the method 
values for deflection the grid joints are tabulated Table 
Since 1,000 was used the equation formation, the effect multiply- 
ing 1,000 was obliterated. All deflection values shown are 
multiples which the and are properties the grid beam. 

example moments and torques may determined after the 
deflections are known. 

Example C—Analysis Uniformly Loaded Square with Clamped Edges.— 
The square slab shown Fig. has been selected for analysis show the 


Special Publication No. 28, Coast and Geodetic Survey, Washington, C., 1915. 
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accuracy the grid method used define the entire elastic surface 
loaded slab. The characteristics this slab are follows: 4,000,000; 
load per unit area. The dashed lines Fig. represent the grid 
beam centerline. Deflection values obtained ten different points the 


ABCDDCBA 


> > 


slab will compared with the exact values obtained the mathematical 
theory 

Nine grid beam centerlines, each two inches apart, are selected the two 
orthogonal directions, order convert the slab into pseudogrid system. 
This slab division gives rise the ten singular joints shown. Under the action 
uniform load, the square slab will deflect symmetrically around its 
south, east-west and diagonal axes. Because this symmetry, the singular 
joints contained within any one octant the grid will completely define the 
elastic surface. 

Since the grid beams both directions have identical lengths and cross 
sections, the resistances bending and twist will constant every joint. 
When the physical characteristics the slab and the constant length and 
cross-sectional properties the grid beams are involved the proper expres- 
the following distribution factors will result: 


Percentage moment 13.49 torque..............(18a) 
Percentage moment 36.51 


unbalanced moment applied any joint the gridwork 
either direction, 36.51% the applied moment resisted two beams 
bending, and two beams torsion resist 13.49% the applied moment. 

The multiple symmetry existent this example permits the simultaneous 
displacement all similarly numbered joints and the calculation the 
auxiliary force systems developed these displacements. example, 
the eight joints designated Fig. will deflected through equal dis- 
tances when the gridwork uniformly loaded. Furthermore, these eight 


Calculation Maximum Deflection, Moment, and Shear for Loaded Rectangular 
Plate with Clamped Wojtaszak, Transactions, American Society Mechanical Engineers, 


Vol. 59, 1937, A-173. 
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joints form one series singular joints. They are singular because their 
displacement, when coupled with prohibition displacement all other 
joints, will develop distinctive moment and torque pattern. 

When all the joints designated are simultaneously displaced, 
shown Fig. 14, symmetrical distribution bending and torsional moments 
will result from the symmetrically located fixed-end moments, Consequently 
the computations can confined the upper left quadrant. 

The numerical process distribution illustrated Fig. for the dis- 
placements joints only. Again, the fixed-end moments are expressed 
multiples and the signs these moments are keeping with 
the usual moment distribution conventions when viewed from the south. 


1045 


Center Line 


Only the east-west, fixed-end moments are considered the grid system be- 
cause the symmetry the structure allows identical force patterns 
developed the north-south direction, when north-south, fixed-end moments 
are distributed. 

After the distribution process complete, final bending moment values 
the grid beam ends are used compute end shear values terms the original 
displacements. The force pattern developed the north-west quadrant 
shown Fig. 16. This force system not complete, however, because 
shows only the effects the east-west, fixed-end moments induced the 
displacements As. Fortunately, the force pattern developed the north- 
south, fixed-end moments can expressed without further calculation. For 
example, the reactions 102, 483, 72, 654, read from 
east west along grid beam shown Fig. 16, will also the north- 
south direction along the orthogonal grid beam The sixth column 
Table 2(a) shows the total reaction values due the simultaneous displace- 
ments. The total reaction values result from algebraic summation 
reaction magnitudes developed both north-south and east-west, fixed-end 
moments. 

Similarly, the remaining groups singular joints are displaced, the 
fixed-end moments distributed, and the auxiliary forces These 


— 


GRIDWORKS AND SLABS 


TABLE EQUATIONS AND SOLUTIONS FOR DEFLECTION, 
UNIFORMLY LOADED SQUARE PLATE 


0.272 0.766 1.395 2.135 5.135 


3.304 3.935 


6.134 7.341 


operations are not included because their repetitious character. All auxiliary 
force values, expressed multiples are shown the left 
side Table 

The ten linear equations defining the elastic surface the slab under 
load are found equating the separate reaction magnitudes each joint 
the externally applied the same joint. Consequently, Table 
not only expression auxiliary force systems but also statement the 
ten linear equations equilibrium. 

The force everywhere equal For convenience, the load 
multiplied 1,000, thus nullifying the effect the thousandfold increase 
moment value. The solutions these equations are shown Table 
which the values displacements are multiples 


Method 


Method 


Grid analogy .00140 .00215 .00255 .00390 .00604 .00720 .00939 .01122 .01343 


this classical slab-problem, the magnitudes deflections over the 
loaded surface, determined the theory elasticity and the grid 
analysis procedure, are found excellent agreement. comparison 
deflections the ten selected positions shown Table The results 
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slab; and the other terms are previously defined. 

Since the length the grid beam this example equal a/9, the 
grid analysis values shown Table can converted follows: 


Along the slab boundary there less grid beam definition the nearest 
grid beam one half grid beam length removed from the edges. result 
the greatest error the grid analysis occurs near the slab boundary. 


CONCLUSIONS 


The analysis gridworks and slabs subjected normal load tedious 
process, for the calculation deflections and the determination their sig- 
nificance should performed many points order define clearly the 
resulting elastic action. Neither the exact mathematical methods nor the 
approximation schemes can circumvent this difficulty. 

The procedure developed this paper should compare favorably with 
any known method for determining slab and gridwork deflections. The method 
fundamentally numerical process that results the solution only one 
series simultaneous equations. When any series linear equations contain 
frequent coefficients zero values for the unknowns, the solution the 
equations relatively simple. This condition may realized the slab 
deflection equations this method performing fewer cycles distribution, 
thus confining the dispersed moments fewer joints. Obviously, more errors 
are incurred when this expediency employed. 

When applied true gridwork, the outlined method results 
solution, for degree indeterminateness avoided. The moment and 
torque distributions express the unknown moments and twists each joint 
terms the unknown displacements that comprise the terms the linear 
equations. believed that the initial threefold indeterminacy assumed 
each joint represents structural condition more realistic than ordinarily 
considered other methods grid analysis. 

While the authors can offer mathematical proof the analogy between 
slab and gridwork, numerous results computed the grid method and checked 
the theory elasticity seem validate the grid analysis procedure. Ob- 
viously, the accuracy any analysis dependent the number grid beams 
used approximate the slab. Using analogous grid consisting six 
north-south beams and six east-west beams and solving directly linear equa- 
tions containing all three unknowns each joint, Béla Enyedi® has shown 
that the computed center moment uniformly loaded square plate agrees 
remarkably well with Henri Marcus’ solution the theory elasticity. 
interesting note that the deflections example computed the grid 
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method, agree more closely with the exact values than those offered Mr. 
Marcus’ finite difference using the same slab division. 

The method moment and torque distribution for determining slab and 
gridwork deflections relatively simple apply wide variety problems, 
heretofore difficult solve. 

The degree and the variation fixity freedom existing the slab 
gridwork boundary does not complicate the solutions for deflections and 
moments, since the distribution factors grid beam terminals are assigned 
reproduce the given boundary condition. When the fixed boundary 
conditions the slab example were imposed the grid beam ends, the 
elastic surface developed the grid system was found strikingly similar 
that computed Mr. Wojtaszak. Fig. illustrates slab irregular 
shape and variable boundary conditions and the possible arrangement 
grid beams defining the various conditions. 


5555555553557 


ABLE BouNDARY CONDITIONS 


The analysis flat slabs that are continuous over column wall supports 
possible this method. shown example the center column support 
actually simplifies the analysis the flat slab prohibiting deflection the 
support joint and consequently reducing the number unknowns. 

Slabs variable cross sections, when divided into analogous grid beams, 
submit solutions the variation distribution factors considered. 
Tentative grid beam layouts for various cases are shown Fig. 18. 

Other extensions the moment and torque distributions are feasible. 
Investigations are being made slabs elastic foundations order 
predict deflections and stresses the slabs. 

general, the moment distribution method for determining deflections 
and moments gridworks and slabs exhibits versatility that not char- 
acteristic other analytical methods. There are some who will reject the 
simplicity the numerical approximation for the elegance the exact mathe- 
matical analysis. Where solutions the theory elasticity are known, this 
attitude seems justified, although thorough evaluation slab action 
any scheme time-consuming process. The real merit the grid analogy 
approximation used this paper has yet explored slabs heretofore 
defying solution. 


Theorie elasticher Gewebe und ihre Anwendung auf die Berechnung biegsamer 
Henri Marcus, Julius Springer, Berlin, 1924, 152. 
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GRIDWORKS AND SLABS 


APPENDIX.—COMPUTATION DEFLECTIONS 


Wojtaszak gives the following expressions for the calculation 
deflection uniformly loaded rectangular plate with clamped edges: 


cos 
cos 
B,, 
cosh 


which the deflection; the load per unit area; the length 
side; and are coordinates plate, origin center; the plate 


n An Bn 


0.0365 0.0365 


The deflection values shown Table were calculated using these expres- 
sions. Deflection influences can determined for any point defined the 
coordinates and expanding the above series. Since the series converges 
very rapidly, only three four terms were deemed necessary. Throughout 
the calculation accuracy five decimal places was used. 
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